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ABSTRACT 



Astronomical mid-IR spectra show two minor PAH features at 5.25 and 5.7 
/im (1905 and 1754 cm~^) that hitherto have been httle studied, but contain 
information about the astronomical PAH population that complements that of 
the major emission bands. Here we report a study involving both laboratory 
and theoretical analysis of the fundamentals of PAH spectroscopy that produce 
features in this region and use these to analyze the astronomical spectra. The 
ISO SWS spectra of fifteen objects showing these PAH features were considered 
for this study, however only four (HD 44179; NGC 7027; Orion Bar, 2 positions) 
have sufficient signal-to-noise between 5 and 6 /im to allow for an in-depth anal- 
ysis. All four astronomical spectra show similar peak positions and profiles. The 
5.25 iJ,m feature is peaked and asymmetric, with a FWHM of about 0.12 ± 0.01 
/im (~40 ± 6.5 cm~^), while the 5.7 /im feature is broader and flatter, with a 
FWHM of about 0.17 ± 0.02 //m (50 ± 5.6 cm-^). Detailed analysis of the lab- 
oratory spectra and quantum chemical calculations show that the astronomical 
5.25 and 5.7 /xm bands are a blend of combination, difference and overtone bands 
primarily involving CH stretching and CH in-plane and CH out-of-plane bending 
fundamental vibrations. The experimental and computational spectra show that, 
of all the hydrogen adjacency classes possible on PAHs, solo and duo hydrogens 
consistently produce prominent bands at the observed positions whereas quar- 
tet hydrogens do not. In all, this study supports the picture that astronomical 
PAHs are large with compact, regular structures. Prom the coupling with pri- 
marily strong CH out-of-plane bending modes one might surmise that the 5.25 
and 5.7 /im bands track the neutral PAH population. However, theory suggests 
the role of charge in these astronomical bands might also be important. 
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Introduction 



Some thirty years of observations, combined with computational and laboratory 
studies, have shown that the mid-IR astronomical emission features, formerly referred to as 
the Unidentified Infrared (UIR) bands, are produced by mixtures of highly vibrationally 
excited Polycyclic Aromatic Hydrocarbons (PAHs) and closely related species. Detected in 



many Ga. 



YanetaL 



actic and extra-galactic objects, including several with significant redshift (e.g. 



20051 ). the astronomical infrared emission features present an important and 



unique probe of astrochemical and astrophysical conditions across the universe. Recen t 



reviews and papers of t 



le observational anc 



Hudgins fc AUamandola 



Tieleni 



Li fc Draine 



2004 



van Dishoeck 



laboratory work (e . 



2004 



2008) and work on theoreti c al models (e.g. 



2001 



Pech et al. 



2002 



Sellgren et al. 



Peeters et al 



2007 



Verstraete et al 



Rapacioli et al. 



2006 



2004 



Smith et al. 



2001 



Draine fc Li 



2007: 



Bakes et al 



2001 



20071 ) can be found 



elsewhere. 

The major features at 3.3, 6.2, '7.7', 8.6 and the complex of bands between 11 an d 



20 um have been s 


tudied 


in great detail (e.g. 


Van Kerckhoven et al. 


2000 




Honv et al. 


2001; 


Peeters et al. 


2002 




van Diedenhoven et al. 


2004 


, and references therein), and the 



fundamental spectroscopic information is now available with which one can analyze the 
strongest astronomical features. However, there are several components of the astronomical 
PAH emission spectra that have been widely overlooked. Many of these contain valuable, 
sometimes subtle, information which is equally important to that revealed by the more 
well-known features. This study focuses on just such features, namely the weak bands that 
fall between 5 and 6 /xm (2000 and 1667 cm~^). 

One of the early predictions of the PAH hypothesis was the expectation of a weak 
emission feature near 5.25 / xm in all objects showing the major PAH bands. Its detection 
in 1989 fi AUamandola et al.i ) was an early confirmation of the PAH hypothesis. Since that 
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time, although evident in many spectra showing the major PAH features, httle has been 
pubhshed on this feature and its companion near 5.7 /^m. 

The weak 5.25 and 5.7 /im PAH features do not correspond to fundamental vibrational 



< Uj), combinations 


+ V 


AUamandola et al. 


1989a) 



frequencies (i/j, z/^, ■ • ■), but are produced by overtones 

difference ivi — i/j)bands of these fundamental vibrations (j AUamandola et al. 
example, the strong CH out-of-plane (CHoop) fundamental bending vibration (t'oop) for solo 
hydrogens produces the well known band at 11.2 yum (893 cm~^). The overtone of this 
vibration, 2 x i^oop, is expected to produce a much weaker feature near 5.6 /im (1786 cm^^) 
and is seen to contribute to the blue side of the broad 5.7 /xm interstellar feature. Likewise, 
z/oop could combine with a CC stretching vibration (t'oop ± '^cc); resulting in other weak 
features. 



Here we present high quality ISO SWS fide Graauw et al. 



19961 ) spectra from four 



astronomical sources that show these features. The 5.25 and 5.7 yum bands are analyzedin 
terms of overtone, combination and difference frequencies using experimental and theoretical 
PAH spectra. This manuscript is structured as follows. The astronomical observations are 
presented in Sect. 121 PAH spectroscopy is described in Sect. [3l astrophysical implications 
are drawn in Sect. [Hand a summary with conclusions is given in Sect. [51 



2. The Astronomical emission features in the 5-6 /im (2000 to 1667 cm ^) 

region 

2.1. Observations 

A sample of high quality spectra from 15 sources, obtained by the Short Wavelength 
Spectrograph (SWS) onboard ESA's Infrared Space Observatory (ISO), were investigated 
for PAH band emission in the 5-6 /im region. Of this set, the spectra from four objects 
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have sufficient signal-to-noise in this wavelength range permit an in-depth analysis. 



The four sources considered here are: HD 44179 ('Red Rectangle'), NGC 7027, and 
two positions on the Orion ionization ridge, H2S1 and D5. HD 44179 is the central binary 
star system of a bipolar planetary nebula, with o ne member being a p ost-AGB star. 



Both stars have a common circumstellar disk (e.g. 



Waelkens et al. 



19961 ). NGC 7027 is a 



compact 



( iLatter et al 



c m) carbon-rich Planetary Nebula with a hot White Dwarf at its center 



2000l ). The Orion Bar spectra probe the Photon Dominated Region forming 



the interface between the bright H II region that is ionized by the Trapezium stars, and the 
Orion Molecular cloud; H2S1 and D5 are two positions within the bar. Table [T] shows the 
journal of observations and the available astrometric data. 



2.1.1. Data Reduction 



The data have been obtained using the AOT 01 full scan mode at various speeds or the 
AOT 06 mode, providing spectra with resolving power of 400 to 1600 (A/AA). The data 
were processed with lA^, the SWS Interactive Analysis package, using calibration ffies and 



procedures equivalent with OL 



be found in 



version 6.0 or 10.1. A detailed account of the reduction can 



Peeters et all (120021 ). Fig. [Dpresents the raw extracted 2.38 - 45.2 /im ISO SWS 



spectra of the four sources at their original resolution. These data have not been defringed. 
The OSIAIJ software package was drawn on to further reduce and analyze the data between 
5 and 6 /im. Further reduction included bad data removal (e.g. glitches), edge truncation, 
splicing and merging of detector bands. Uncertainties were estimated by comparing the up 
and down scans, yielding a measure for the noise. For HD 44179 a difference in absolute 



''OSIA is a joint development of the SWS consortium. Contributing institutes are SRON, 
MPE, KUL and the ESA Astrophysics Division 
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flux levels, not in shape, between the up and down scan was found. The up scan was scaled 
such that the mean of the up scan would match the mean of the down scan to compensate 
this. This can be done because our interest lies only in the uncertainty in relative flux 
level, since those trace the reality of the shape of a feature. When required, the 5-6 /im 
spectra were binned to a constant resolution of 800 to facilitate one-to-one comparison. For 
Orion Bar position D5 this results in oversampling since the original data were obtained at 
a resolution of about 400. Table [2] presents a subset of the data reduction parameters and 
Fig. [2] presents the fully reduced 5-6 /im ISO SWS spectra of the four sources at their 
original resolution. 



2.2. Analysis 



The spectra display a wealth of structure, including the familiar PAH bands at 3.3, 
6.2, '7.7', 8.6, 11.2 and 12.7 /im with their underlying plateaus and sub-features. For 
clarity, several atomic emission lines, have been masked out in Fig. [H Detailed analysis 



of the major 



Peeters et al 



eatures in the spectra have been p ubhshed elsewhere (e.g. 



2002 



van Diedenhoven et al. 



Hony et al. 



2001 



2004J ). Here we focus on the PAH emission 



between 5-6 /im. As Fig. [T] illustrates, the bands in this region are weak compared to the 
major features. 

To isolate the PAH band profiles a continuum is established for each object. The 
continua used apply only to the 5-6 /im region. Spline continua were constructed using 
connecting points between 5-6 /im, straddling the 5.25 and 5.7 fim features. The adopted 
continua are shown in Fig. [2] and the connected wavelengths are listed in Table O 

Atomic emission lines can also contribute to the emission in this region. For some of 
these objects narrow lines at 5.128 (1950 cm'^; H I 6 - 10), 5.510 (1815 cm'^; H2 
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Table 1. Journal of observations and available astrometric data. 



Source 


a (2000) 


<5 (2000) 


TDT 


Obs. mode 


Sp. 


Go 


Object Type 


Ref. 




[h m s] 


[° ' "] 




(speed) 


Type 


[Habing] 






HD 44179 


06 19 58.20 


-10 38 15.22 


70201801 


01 (4) 


B8V 


5 ■ 10^ 


Post- AGE star 


1 


NGC 7027 


21 07 01.70 


+42 14 09.10 


55800537 


06 


2 10^ K 


2 ■ 10^ 


PN 


1 


Orion Bar H2S1 


05 35 20.31 


-05 25 19.99 


69501806 


01 (4) 


06 


7 - 10^ 


H II region 


2 


Orion Bar D5 


05 35 19.81 


-05 25 09.98 


83101507 


01 (2) 


06 


5-10'' 


H II region 


3 



References. — (1) 



Beintema et al 



1 19961) : (2) 



Verstraete et al 



i 200lh : (3) 



Peeters et al 



12003) 



Table 2. Data reduction parameters and connected wavelengths. 



Source 


OLP 


Resolution 


Continuum points 


Integration bounds 












HD 44179 


10.1 


-800 


5.017;5.494;5.836;5.922 


5.160;5.355;5.580;5.830 


NGC 7027 


6.0/10.11 


-1600/1000 


5.025;5.079;5.832;5.938 


5.150;5.375;5.570;5.805 


Orion Bar H2S1 


10.1 


-800 


5.102;5.415;5.825;5.969 


5.170;5.374;5.590;5.820 


Orion Bar D5 


10.1 


-400 


5.057;5.382;5.472;5.943 


5.150;5.380;5.560;5.845 



References. — (1) The AOT 01 and AOT 06 data have been reduced with OLP version 6.0 and 10.1 
respectively. The AOT 01 data is used in Fig. [T] and the better quahty AOT 06 data is presented in Fig. [2] 
and used onwards. 



WQvenumber [cm ^1 

JOOO 2000 1500 1100 1000 900 800 700 600 




4 6 8 10 12 14 16 

wavelength [/xm] 



Fig. 1. — Extracted mid-IR spectra from 3 - 14.5 /xm of the four sources considered here: HD 
44179, NGC 7027 and the two positions along the Orion Bar. Spectra arc shown at original 
resolution. For clarity the prominent atomic emission lines, e.g. at 5.61 ([Mg V]), 10.51 
([S IV]) and 12.81 ([Ne II]) //m, have been masked out. The arrows indicate the positions of 
the subtle 5.25 and 5.7 /im features. 
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wavenumber [cm ^] 

2000' 950' 900 1 850 1 800 1 750 170 




5.0 5.2 5.4 5.6 5.8 6.0 
wavelength [fJ^m] 



wQvenumber [cm ^] 




... I ... I ... I ... I ... 1 
5.0 5.2 5.4 5.6 5.8 6.0 
wavelength [fim] 



Fig. 2. — The fully reduced 5-6 /xm spectra of the four objects considered. The are presented 
at original resolution. Also shown are the adopted continua. 
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S(7)) and 5.61 /iin (1783 cm~^; [Mg V]) are present. Special care was taken to remove the 
5.61 fim [Mg V] line in the spectrum of NGC 7027 where it is blended with the 5.75 fim 
PAH feature. The resulting, continuum subtracted spectra are presented at equal resolution 
in Fig. [3l Two broad bands, centered roughly at 5.25 and 5.7 /im, are clearly present in 
all objects and are quite similar within the uncertainty (systematic + noise). The 5.25 yum 
band is asymmetric and sets off near 5.15 /im (1942 cm~^). For both positions in the Orion 
Bar it also seems to start off with a satellite feature, centered at about 5.2 fim (1923 cm~^) 
and similar to those seen at 6.0 and 11.0 fim for the 6.2 and 11.2 fim bands, respectively. 
After a steep rise the profile peaks at 5.25 fim (1905 cm~^) and slowly declines to reach 
zero intensity near 5.4 /xm (1852 cm"^). This red wing resembles that seen for the 3.4, 6.2 



and 11.2 fim features (IRoche et al. 



19961 ). Fitting the profiles at original resolution with a 



single Gaussian results in a FWHM of about ~0.12 ± 0.01 fim (~42 ± 6.5 cm ^. 

The 5.7 fim feature is broader, and starts its rise near 5.55 fim (1802 cm~^), just 
beyond the 5.51 fim (1815 cm^^) H2 ^ S(7) line. The profile is unusual from about 5.65 
fim (1770 cm^^) to about 5.75 fim (1739 cm^^), where the feature seems flat-topped and 
then drops to the continuum at 5.85 fim (1709 cm^^). Some spectral structure is visible in 
on the fiat-top, giving the impression of a blended double-peaked feature. Fitting the profile 
with a single Gaussian results in a FWHM of about 0.17 ± 0.02 fim (~51 ± 5.6 cm^^). 

After emission line removal the bands are integrated in the A — Fa domain, using a 
multiple Simpsons' rule, yielding band strengths in W-m~^. Uncertainties are calculated 
from the data. Table [3] summarizes these results and compares them with the band 
strengths of the dominant PAH bands. 
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wavenumber [cm"^] 
2000 1950 1900 1850 1800 1750 

10 I \ \ \ \ \ — 

HD 44179 



1700 



8 




NGC 7027 



6 




Orion Bar H2S1 




Orion Bar D5 



Galaxies 








5.0 



5.2 5.4 5.6 5.8 

wavelength [micron] 



6.0 



Fiff. 8.— The 5.25 and 5.7 urn PAH features in HD 44179. NGC 7027 and the Orion Bar: 
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3. PAH Spectroscopy in the 5 to 6 fim (2000 to 1667 cm^) Region 

3.1. Experimental Studies 

The matrix isolation infrared s pectroscopy techniques employe d in this study haye 



been described in detail previously (IHudgins fc AUamandola 



1995a 



a; 



Hudgins &: Sandford 



19981 ) and will be summarized here only briefly. 



Matrix isolated PAH samples were prepared by vapor co-deposition of the species 
of interest with an overabundance of argon onto a 14K Csl window suspended in a 
high- vacuum chamber (p < 10~^ Torr). The samples were vaporized from heated Pyrex 
tubes while argon was admitted through an adjacent length of copper tubing, cooled by 
liquid nitrogen. The deposition temperatures, which are dependent on size and structure as 
well as the sources of the i ndividual PAHs, are given in the references summarized in Table 
1 of 



Mattioda et al. 



(120051 ). 

Spectra from 1.7 - 20 /im (6000 to 500 cm~^) were measured on either a Nicolet 740 or 
a Digi-Lab Excalibur FTS 4000 FTIR spectrometer using a KBr beam splitter and a liquid 
nitrogen cooled MCT detector. Each spectrum represents a co-addition of between 500 and 
1024 scans at a resolution of 0.5 cm^^. The number of scans was chosen to optimize both 
the signal-to-noise as well as time requirements for each experiment. 

Figure m shows the 5 to 16.5 spectra of several matrix isolated neutral PAHs. These 
spectra are typical of all spectra in the Ames collection of neutral PAH IR spectra and 
serve to illustrate that all PAHs have weak bands in the 5 to 6 /im region. For presentation 
purposes only, some of the data has been baseline corrected using the Digi-Lab data analysis 
software package. No further data reduction was necessary. All numerical values were 
obtained from the original (unaltered) data. 
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wavelength (i^m) 

5 5.2 5.4 5.6 5.8 9 10 11 12 13 14 15 16 




1950 1800 1200 1100 1000 900 800 700 600 

wavenumber (cm"^) 



Fig. 4. — Matrix isolated spectra for the CHoop region (10 - 15 //m) as well as the overtone 
and combination band region (5-6 ^m) for various neutral PAHs containing solo, duo, trio 
and quartet hydrogens. The number of adjacent hydrogen atoms per ring are indicated in 
the structure. 
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3.2. Theoretical 



2003 



Hudgins et al. 



Previous work (e.g. iBauschlicher fc LanghofllllQQTl : iBauschlicherl l2002l : iMattioda et al 



2005 



Bauschlicher et al 



20081 ) h ave shown t 



theory (DFT) calculations performed using the hybrid (IBecke 



19941 ) approach in conjunction with the 4-3 IG basis set (IFrisch et al.lll984l ) yield scaled 



lat densit y functional 



1993h B3LYP fIStephens et al. 



harmonic frequencies that are in excellent agreement with results obtained from matrix 
isolation experiments. In this work, we expand our previous work by explicitly considering 
anharmonic effects (fundamentals, overtones, combination, and difference bands), which 
are computed by doing numerical differentiation along the normal modes 



Barone 



20051) 



after the geometries are fully optimized and the harmonic frequencies are computed using 



analytic second deri vatives. All o : 



03 progra m system (IFrisch et al. 



Molekel 



(jPortmann fc Liithi 



the DFT calculations are performed using the Gaussian 



2000). 



20041 ). The vibrational modes are illustrated using 



As discussed previo usly, the DFT results differ sys tematically from experiment and 



hence need to be scaled ( IBauschlicher fc LanghofI 



19971 ). The scaling factor depends on the 



level of theory (the correlation treatment and basis set), and it is common that CH stretches 
will have a different scaling factor from the other modes because of the large anharmonic 
effects for hydrogen stretching motions. In this work we use the small 4-31 G set that has 
been used successfully in previously and t he significantly lar ger and more precise, but also 



more time consuming, 6-311G** basis set (IFrisch et al. 



19841 ). The larger 6-3 IIG** basis set 



was used for naphthalene (CioHg) because of its small size. The results from the large basis 
set allow us to calibrate the smaller ba sis set that is required for the larger molecules. We 



use the matrix isolation experiments of 



Hudgins et al. 



(119941 ) for the calibration of neutral 



^MOLEKEL 4.3, P. Fliikiger, H.P. Liithi, S. Portmann, J. Weber, Swiss Center for Sci- 
entific Computing, Manno (Switzerland), 2000-2002 
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naphthalene. While gas-phase data exists, we use the matrix isolation experiments to allow 
a direct comparison between our experim ents and our computed value s . In a ddition, it has 



been shown by e.g. 



Oomens et al. 



fl2000[ ) and 



AUamandola &: Hudgind (120031 ) that ther is a 



very good agreement between the gas phase and matrix isolation results. 

Comparison of the computed fundamental frequencies of neutral naphthalene, using 
the 6-311G** basis set, shows that a scale factor of 0.998 brings the computed non-CH 
stretching modes (i.e. all the bands excluding the CH stretching modes) into good 
agreement with experiment. This comparison shows an average absolute error of 2.6 cm^^ 
and a maximum absolute error of 9 cm~^ between experiment and theory. A scale factor of 
1.005 brings the computed CH stretching modes into good agreement (a maximum error of 
7 c m^^) with experim ent. This accordance with experiment is consistent with recent work 



by 



Cane et al 



(120071 ). who find good agreement between theory and experiment using a 
hybrid functional and a large basis set (B971/TZ2P). For the smaller 4-31 G basis set we 
find a scale factor of 0.976 for the non-CH stretching fundamentals and 1.007 for the CH 
stretching fundamentals. The maximum absolute error in the non-CH stretching modes 
is 19 cm"'^, a value greater than that when using the larger basis set, but still sufficiently 
small to allow a qualitative understanding using the small basis set. It is interesting to note 
that two scale factors are needed for the 4-31G fundamentals while only one scale factor is 
needed for the 4-31 G harmonics. 

Since the scale factors for the 6-31 IG** basis set are very close to 1, we do not use 
any scaling for the overtone, combination, or difference bands. For the 4-31G basis set, 
the overtone and combination bands are scaled by 0.976, the scale factor for the non-CH 
stretching modes, because the CH stretches are not involved in the overtone or combination 
bands of interest. For the difference bands, the non-CH stretching fundamentals are scaled 
by 0.976 while the CH stretching fundamentals are scaled by 1.007. While the level of 
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theory used is expected to yield accurate band positions for the overtones, combination, 
and difference bands, the current version of Gaussian does not determine intensities for 
these bands. However, as discussed in Sec. 13. 3. 2^ selection rules allow the exclusion of 
some combination and difference bands as possible contributors to the observed features. 
In addition to excluding some bands, we speculate on which allowed bands will be the 
strongest. In the harmonic approximation, the second derivative of the dipole moment 
is responsible for the intensity of the first overtone, but in the real world, anharmonicity 
means that the overtones will have some contributions from the first derivative. This implies 
that the strongest overtone, combination, and difference bands will commonly involve the 
strongest fundamentals. 

Figure [5] compares the 5 to 19 fim (2000 to 530 cm^^) spectrum of neutral naphthalene 



computed using density functional theory at the B3LYP l evel t o the spectrum of matrix 



isolated neutral naphthalene (CioHg) from 



Hudgins et al. 



(119941 ). This figure demonstrates 



the excellent agreement between the theoretically computed and experimentally measured 
mid-IR spectrum for fundamentals and the reasonable agreement in the positions of the 
calculated overtone and combination bands. Note that difference bands only arise from 
emission processes and are thus not found in the experimental spectra. 



3.3. The Origin of the PAH Bands in the 5 to 6 /xm (2000 to 1667 cm ^) range 



The experimental spectra presented in Fig. H] illustrate what has long been known, 
PAHs show weak spectral structure in the 5 to 6 fim (2000 to 1667 cm^^) region. This 
structure has been attributed to ove rtones and combin ation bands based on spectroscopic 



studies of small aromatic molecules (lYoung et al. 



19511 ). This work is summarized below. 



followed by a detailed analysis of the PAH vibrations that produce bands in the 5 to 6 /xm 
region for large PAHs, species comparable in size to the astronomical PAHs. 
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J I L 
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Li 
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J I L 



J I l_ 



J I L 



J I l_ 



2000 1800 1600 1400 1200 1000 

wavenumber (cm"^) 
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Fig. 5. — Comparison between the B3LYP-level calculated and experimentally obtained 
spectra of naphthalene from 5 to 19 fim, showing their outstanding agreement for funda- 
mentals. The agreement for the overtone and combination bands is reasonable. *Note that 
since the theoretical results for the overtone and combination bands do not calculate inten- 
sities. Therefore, the bands were derived from fundamentals with an intensity >10 km/mole 
and given the same intensity to allow for comparison with the experimental frequencies. 
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Young et al. 



(I1951I ) present some examples of the distinctive band patterns in the 5 to 
6 fim region produced by different substituted benzenes. The different substitution patterns 
force different H adjacencies and this has been found to determine the characteristic 
band patterns in this region. This behavior strongly suggests that fundamental vibrations 
associated with hydrogen are the primary carriers of this spectroscopic structure, but 
these different substitution patterns could also induce IR activity in transitions involving 
CC vibrations. The central role of hydrogen in producing these features is confirmed 
when one compares the 5 to 6 /im region of the spectra of normal PAHs (CxHy) with the 
same region in the spectra of fully deuterated PAHs (PADs; CxDy). The absence of any 
activity in the 5 to 6 /im region of the deuterated PAH spectra shows that hydrogen is 
essential in producing these features. The band patterns in the 5 to 6 fim region reflect 
hydrogen adjacency because, as discussed below, they involve strong CH vibrations whose 
fundamental frequencies are determined by their adjacency class. 



3.3.1. Experimental 



The 5 to 6 /im region of the astronomical emission spectrum was initially thought 
to be dominate d by simple combinations and overtones of the fundamental CHoop modes 



(IBellamy 



19601 ). Thus our analysis begins by separating out the contributions of the CHoop 
modes. The CHoop region, which falls between 10.5 and 15 fim, includes contributions of 
the solo (non-adjacent hydrogens), duo (doubly adjacent hydrogens), trio (triply adjacent 
hydr ogens) and quartet (quadruply adjacent hydrogens) hydrogens present on a PAH 



[e.g. 



Hudgins fc AUamandola 



19991 ). Each of these type of vibrations falls within specific 



wavelength regions. Given the segregation of the CHoop vibrations in PAH spectra one 
would expect combination modes arising from these fundamental vibrations to be similarly 
grouped (i.e. quartets furthest to the red and solos furthest to the blue), if the various 
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CHoop modes combined with the same vibrational modes. To better understand this region, 
we have analyzed the spectra of over 35 neutral PAHs of various sizes and structures and 
have found that this is not necessarily the case. While several features between 5 and 6 /xm 
appear to be the result of solo, duo, trio and quartet hydrogens, their order does not follow 
that observed in the CHoop region. 

The CHoop spectra of several matrix isolated PAHs are shown in conjunction with 
their respective overtone and combination regions in Fig. HI As one can easily see, the 
CHoop region varies considerably with respect to the various PAH structures, reflecting the 
variability in the solo, duo trio and quartet hydrogens from one PAH to the next. Likewise 
there is considerable variability in the overtone and combination region (5.1 to 5.8 fim 
region). Close inspection of Fig. H] reveals that molecules b, c, and e contain a large set 
of bands spanning the full wavelength region, while molecules a and d have a pronounced 
break in the overtone and combination bands between the 5.25 and 5.7 fim regions. The 
major structural difference in these PAHs is the presence of quartet hydrogens in molecules 
b, c, and e. The presence of quartet hydrogens in a PAH molecule produces a complex of 
bands spanning the entire 5 and 6 fim region, making them unsuitable as responsible for 
the astronomical 5.25 and 5.7 fim features. 

Ignoring, for the moment, the fundamental vibrational modes as carrier of the overtone 
and combination bands, several additional interesting trends appear in the 5 to 6 fim 
region which can be linked to the PAH structure. The PAH spectra presented in Fig. [6] 
will assist in identifying these trends. The positions of the astronomical 5.25 and 5.7 fim 
have been indicated by the shaded areas, with on b oth sides an uncer tainty of ~10 cm~^ 



19951 ). The first trend to 



to account for cold absorption versus hot emission (iJoblin et al.l 
note is that all the matrix-isolated PAHs in Fig. [6l along with all the PAHs investigated 
by this lab to-date, exhibit bands in the 5.25 and 5.7 fim regions. The quartet containing 
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PAHs, in addition to exhibiting bands in undesirable positions, also tend to have bands to 
the blue and red of the positions occupied by the 5.25 and 5.7 fim astronomical features. 
For instance, in Fig. the bands extend further to the blue of the 5.25 fim area than the 
bands visible in the PAHs shown in panels e, f, g and h of Fig. [6l which contain only solo, 
duo and trio hydrogens. Furthermore, quartets that are part of a 'bay' region in PAHs, 
such as found in panels c or d of Fig. El tend to produce bands at wavelengths shorter than 
those of less hindered quartets, such as shown in panels a and b of Fig. [61 Typically the 
bay region quartets result in a band around 5.128 fim (1950 cm~^) versus the 5.14 fim (1940 
cm~^) position of the less hindered quartet. Thus PAHs containing quartet hydrogens will 
exhibit bands outside the area occupied by the 5.25 and 5.7 fim astronomical features. 

Additional analysis reveals that PAHs with duo hydrogens tend to exhibit a band 
around 5.19 /im (1925 cm~^), see panel e in Fig. O The presence of trio hydrogens, in 
addition to duos, tends to produce another feature to the blue of the duo band around 
~5.18 yum (1930 cm^^). This is clearly illustrated in the spectra shown in panels e, f, g 
and h of Fig. El where the trios, labeled as 3, in panel f, g and h of Fig. El produce a 
band further to the blue of the duo band, labeled as 2, in panel e of Fig. El However, 
in the absence of duo hydrogens, the trios produce a single feature near the position of 
the lone duo hydrogens. Finally the presence of hydrogen trios also tends to produce a 
feature just to the red of the 5.25 fim band area (~5.39 fim or 1855 cm^^). This feature is 
readily observed (labeled with a 3) in the matrix isolated spectrum shown in panel f, as 
well as in the spectra shown in panels g and h of Fig. El Although all of the PAHs with 
trio hydrogens exhibited this feature, the intensity of this band depends on the specific 
PAH structure. Thus, PAHs containing only duo and trio hydrogens match the 5.25 /im 
astronomical feature well, but PAHs containing trio hydrogens match only marginally due 
the presence of a band around 5.39 fim. This implies a limited concentration of PAHs 
containing trio hydrogens in astronomical objects. 
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Fig. 6. — Matrix isolated spectra of (a) naphthalene, (b) pentacene, (c) chrysene, (d) 
large, uncondensed PAH containing solo, due, trio and quartet hydrogens, (e) coronene. 
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The matrix isolated PAH absorption spectra in the 5.7 fim astronomical band region, 
running from ~5.5 to 5.8 fim, exhibit considerably more variability than the 5.25 fim region. 
Just as in the 5.25 fim region, every PAH investigated shows bands in the 5.7 /im region. 
However there seems to be a greater degree of segregation concerning the absorptions in 
this region. While the presence of quartets hydrogens produce bands in the 5.7 fim region, 
their presence also tends to produce bands to the blue and red of the region. This is 
exemplified in the spectra shown in panels a, b and c of Fig. [61 Although PAHs containing 
only duo and trio hydrogens match the 5.7 fim feature very well, their bands tend to cluster 
around 5.6 /im (1780 cm^^), whereas the presence of solos extends this structure towards 
5.8 fim (1724 cm^^). This is beautifully illustrated in the PAH spectra shown in panels f, 
g and h of Fig. [61 The PAHs producing the spectra in panels f and g of Fig. [61 contain 
only doubly and triply adjacent hydrogens, while the top spectrum is of a similar PAH 
but containing solo hydrogens. As is easily seen in Fig. [6l although PAHs containing duo 
and trio hydrogens can contribute to the 5.7 /im astronomical feature, solo hydrogens are 
required to match the full 5.7 /xm feature. 

In addition to possessing bands to the blue of the 5.25 /im feature and bands to the blue 
and red of the 5.7 fim feature, PAHs containing quartet hydrogens exhibit an additional 
trait which prohibit them from effectively matching the 5.25 and 5.7 fim astronomical 
regions. This concerns the region between the 5.25 and 5.7 /im astronomical features. As 
shown in panel a, c and d of Fig. [6l PAHs containing quartets exhibit more intense activity 
in the area between the 5.25 and 5.7 /im features than PAHs containing only solos, duos 
and trios. In fact, as all the spectra in Figs. [H and [6l reveal, PAHs containing quartets seem 
to exhibit a continuum of bands throughout the region that includes the 5.25 and 5.7 /im 
features. This is exemplified in the PAH spectrum shown in panel d of Fig. [6l Although 
it is a fairly large PAH, its open structure contains solo, duo, trio and quartet hydrogens. 
Likewise its spectrum exhibits a continuum of bands throughout the region containing the 
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5.25 and 5.7 //m features. 

Two final notes should be made regarding PAH structure and the overtone and 
combination region. First, not all of the experimentally measured features in this region can 
be traced back to the strongest solo, duo, trio and quartet hydrogen bands. This precludes 
a one-to-one assignment with the theoretical results for all the bands presented in the 
next section. This region also contains bands around 1903 cm^^ (5.25 iim; ±5 cm~^) and 
1915 cm~^ (5.22 //m; ±5 cm~^) which appear in the PAH spectra collected, regardless 
of structure. In a few instances either the 1903 or 1915 cm~^ band will appear by itself, 
but in most instances both bands are present. Due to the appearance of at least one of 
these bands in all of the PAH spectra, we were unable to associate them with a particular 
structural feature. However, the absence of any bands between 5-6 fim in fully deuterated 
PAH spectra indicate that hydrogen motion is involved, obviously through CH in-plane 
(CHip) modes. Second, while difference bands only arise from emission processes and are 
thus not found in the absorption spectra presented here, they will likely contribute to the 
astronomically observed features. 

In summary, the experimental evidence suggests that the astronomical 5.25 and 5.7 
yUm features can be accommodated by the overtone and combination bands resulting from 
PAHs containing trio, duo and solo hydrogens. The presence of a band around 5.39 /im in 
PAHs containing trio hydrogens implies a hmit of the concentration of this particular PAH 
structure in the astronomical PAH population. The presence of quartet hydrogens produces 
undesirable features in the spectrum. The 5.7 /^m region exhibits a greater variability in 
structure than the 5.25 //m. This greater variability is the result of PAHs with only duo 
and trio hydrogens producing absorption features which cluster on the blue side of the 5.7 
Hm region, whereas the presence of solo hydrogens is necessary to produce features towards 
the red side of this region. Such variability could be utilized as a tracer of PAH structures 
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in astronomical environments. 

3.3.2. Theoretical Studies 

Computational spectra permit an extensive analysis of the nature of the transitions that 
produce features in the 5-6 fim region, including difference bands which are inaccessible 
through experimental absorption spectroscopy. Figure [5] compares the spectrum of neutral 
naphthalene computed using density functional theory at the B3LYP level to the spectrum 
of matrix isolated neutral naphthalene (CioHg). This figure illustrates two important 
points. First it demonstrates the excellent agreement between the experimentally measured 
and theoretically computed fundamental transitions that produce the mid-IR spectrum of 
naphthalene. Second, this level of theory does not predict any fundamental PAH vibrations 
between 5 and 6 fim, adding further support of their assignment to combination and 
overtone types of bands. 

We now consider neutral naphthalene using the 6-3 IIG** basis set. The results are 
summarized in the top two sections of Table |H There are no overtones between 5-6 fim 
that are derived from fundamentals with intensity greater than 10 km/mole. Therefore 
we conclude that overtones are not responsible for the bands between 5 and 6 fim in the 
spectrum of naphthalene. However, naphthalene is a special case since it has only quartet 
hydrogens and, as shown in Table O for coronene, overtones can produce important bands 
in this region. Now, consider combination and difference bands, not overtones. Because of 
the high symmetry in naphthalene, dip ole- allowed combination and difference bands can 
only arise from a forbidden mode interacting with an allowed mode. In addition to this 
rigorous selection rule, we presume that the stronger bands will result from combination 
and differences in which the allowed band is intense (Sect. 13.21) . The combination and 
difference bands of the most interest in the 5-6 fim range are also summarized in Table |H 
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All of the computed combination bands in neutral naphthalene involve the strong CHoop 
bending mode (Bsu) at 12.703 fim. The vibration producing this strong band couples with 
forbidden in-plane and out-of-plane CH bending vibrations to produce the combination 
bands. The fundamental modes are visualized in Fig. [71 Hence, since the only strong bands 
in neutral naphthalene are the CHoop bend and the CH stretch, and the higher frequency 
CH stretch cannot contribute to the combination bands in the 5-6 /im region, the only 
bands in the correct region that can couple with the strong CHoop band to produce 
features in the correct region are other CH bends. Given that difference bands will involve 
a higher frequency mode and that the only strong, allowed higher frequency fundamental 
is the allowed CH stretch, all difference modes falling between 5 and 6 fim must involve 
coupling of the allowed CH stretch with forbidden CHip bending vibrations. 
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Table 3. The integrated band strengths (10 W m ^) of the PAH features between 3 - 
13 /im of HD 44179, NGC 7027 and the Orion Bar: positions H2S1 and D5. Uncertainties 

are given in parentheses in percentages. 



Source 




5.25 


5.75 


6.2I 


7.6l-» 


7 gl,a 


8.6^ 


11.2^ 


12.7 




[nm] 


[Mm] 


[nm] 


[Mm] 


[Mm] 


[Mm] 


[Mm] 


[Mm] 


[Mm] 


HD 44179 


69 (1.4%) 


11 (1.5%) 


19 (1.3%) 


280 (10%) 


121 (10%) 


290 (10%) 


97 (10%) 


113 (10%) 


26 (5.4%)^ 


NGC 7027 


24 (1.7%) 


9.1 (3.3%) 


13 (7.0%) 


95 (10%) 


40 (10%) 


63 (10%) 


28 (10%) 


130 (10%) 


36 (6.9%)^ 


Orion Bar H2S1 


8.2 (1.6%) 


2.9 (1.6%) 


4.4 (1.7%) 


34 (10%) 


32 (10%) 


22 (10%) 


13 (10%) 


29 (10%) 


13 (10%)^ 


Orion Bar D5 


11 (0.55%) 


3.7 (3.2%) 


6.1 (2.3%) 


12 (10%) 


40 (10%) 


36 (10%) 


24 (10%) 


40 (10%) 


17 (7.6%) 


''7.6 and 7.8 /Jm 


bands liavc been 


separated 


using Gaussians 














References. — fl)|peeters et al.l ( 


I2OO2I1: (2)1 


van Diedenhoven 


et al.l J200^ 


a): (3)lHonv et al.l I2OO1I) 
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Table 4. The computed combination and difference bands for naphthalene. Band 
positions determined using the 6-31 IG** basis set are compared to those determined using 
the 4-31G basis set. The forbidden- and allowed transitions that produce each band are 
given. Visualization of the fundamental vibrations producing the combination and 
difference bands indicated here in columns 4 and 7 are shown in Fig. [7] and [HI 

Forbidden transition Allowed transition 

Wavelength Wavelength Symmetry Fig. Wavelength Symmetry Fig. Intensity 
[/im]/[cm~-'-] [/xm] [/xm] [Km mole""^] 



naphthalene 6-3 IIG** neutral 



Combination 








5.989 / 1670 


11.330 


Bag 


b 


5.770 / 1733 


10.570 


Big 


c 


5.660 / 1767 


10.194 


Bag 


d 


5.496 / 1820 


9.680 


Ag 


e 


5.104 / 1959 


8.537 


Ag 


f 


Difference 








5.977 / 1673 


7.324 


Ag 


h 


5.583 / 1791 


8.017 


Bsg 


i 


5.356 / 1867 


8.537 


Ag 


f 


5.309 / 1884 


8.660 


Bsg 


j 


5.894 / 1679 


7.324 


Ag 


h 


5.510 / 1815 


8.017 


Bsg 


i 


5.289 / 1891 


8.537 


Ag 


f 


5.243 / 1907 


8.660 


Bag 


j 



12.703 



3.291 



3.266 



B 



3u 



Biu 



B 



2u 



114.1 



64.1 



49.1 
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Table 4 — Continued 



Forbidden transition Allowed transition 

Wavelength Wavelength Symmetry Fig. Wavelength Symmetry Fig. Intensity 
[//m]/[cm~-^] [//m] [//m] [Km mole~^] 



naphthalene 6-31 IG** cation 



Combination 



5.948 / 1681 


21.365 


Bag 


D 


8.239 


B2u 


A 


218.8 


5.896 / 1696 


10.751 


B2g 


E 


13.047 


Bsu 


B 


105.7 


5.804 / 1723 


19.626 


Ag 


F 


8.239 


B2u 


A 


218.8 


5.779 / 1730 


10.366 


Big 


G 


13.047 


Bsu 


B 


105.7 


5.610 / 1783 


9.826 


B2g 


H 


13.047 


Bsu 


B 


105.7 


5.511 / 1815 


9.541 


Ag 


I 


13.047 


Bsu 


B 


105.7 


5.129 / 1950 


23.198 


B2g 


J 


6.569 


Blu 


C 


87.41 


5.116 / 1955 


13.499 


Big 


K 


8.239 


B2u 


A 


218.8 


5.110 / 1957 


8.410 


Ag 


L 


13.047 


Bsu 


B 


105.7 


5.063 / 1975 


13.131 


Ag 


M 


8.239 


B2u 


A 


218.8 


5.027 / 1989 


21.365 


Bag 


D 


6.569 


Blu 


C 


87.4 


Difference 
















5.415 / 1847 


3.267 


Ag 


N 










5.410 / 1848 


3.266 


Bag 





8.239 


B2u 


A 


218.8 
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Table 4 — Continued 



Forbidden transition Allowed transition 

Wavelength Wavelength Symmetry Fig. Wavelength Symmetry Fig. Intensity 
[/xm]/[cm~^] [/xm] [/^m] [Km mole~^] 



5.371 / 1862 


3.252 


Bag 


P 




naphthalene 4-31G neutral 


Combination 














5.977 / 1673 


11.296 


B2g 










5.776 / 1731 


10.600 


Big 










5.637 / 1774 


10.134 


B2g 




12.677 


Bau 


■■• 111.2 


5.521 / 1811 


9.768 












5.078 / 1969 


8.479 


Ag 










Difference 














5.796 / 1725 


7.386 


Ag 










5.484 / 1823 


7.965 


Bag 




3.248 


Biu 


■■• 78.1 


5.264 / 1900 


8.479 


Ag 










5.212 / 1919 


8.618 


Bag 










5.738 / 1743 


7.386 


Ag 










5.431 / 1841 


7.965 


Bag 




3.229 


B2u 


■■• 71.4 
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Table 4 — Continued 



Forbidden transition Allowed transition 

Wavelength Wavelength Symmetry Fig. Wavelength Symmetry Fig. Intensity 
[/im]/[cm~^] [//m] [/xm] [Km mole~^] 



5.216 / 1917 8.479 Ag 
5.165 / 1936 8.618 B,,^ 
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Table 5. The computed overtone, combination and difference bands for neutral coronene 
computed using the 4-31G basis set. The forbidden- and allowed transitions producing 

each band are given. 

Forbidden transition Allowed transition 

Wavelength Wavelength Symmetry Wavelength Symmetry Intensity 
[/im]/[cm~^] [/im] [/xm] [Km mole~^] 



Overtone 



5.770 / 1733 
Combination 



5.857 / 1707 


11.890 


Big 


5.723 / 1747 


11.349 


Bag 


5.640 / 1773 


11.041 


Bag 


5.475 / 1826 


10.409 


Big 


5.470 / 1828 


10.394 


Bag 


5.432 / 1841 


10.252 


Bag 


5.383 / 1858 


10.091 


Ag 


5.281 / 1894 


9.739 


Ag 


Difference 






5.917 / 1690 


7.206 


Bag 


5.916 / 1690 


7.207 


Ag 


5.732 / 1745 


7.501 


Ag 


5.460 / 1832 


8.024 


Bag 


5.421 / 1845 


8.110 


Ag 


5.399 / 1852 


8.159 


Bag 



11.538 Bau 175.6 



11.538 Bsu 175.6 



3.249 Bi„ 140.3 
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Table 5 — Continued 



Forbidden transition Allowed transition 

Wavelength Wavelength Symmetry Wavelength Symmetry Intensity 
[/xm]/[cm~^] [//m] [//m] [Km mole~^] 



o.oyo 


/ 
/ 


ioOO 


Q 
O 


1 P.P. 

iOO 


A 


5.253 


/ 


1904 


8 


517 


Ag 


5.252 


/ 


1904 


8 


519 


Bsg 


5.989 


/ 


1670 


7 


024 


Bsg 


5.987 


/ 


1670 


7 


026 


Ag 


5.863 


/ 


1706 


7 


206 


Bsg 


5.862 


/ 


1706 


7 


207 


Ag 


5.681 


/ 


1760 


7 


501 


Ag 


5.414 


/ 


1847 


8 


024 


Bsg 


5.375 


/ 


1860 


8 


110 


Ag 


5.354 


/ 


1868 


8 


159 


Bsg 


5.351 


/ 


1869 


8 


166 


Ag 


5.210 


/ 


1919 


8 


517 


Ag 


5.209 


/ 


1920 


8 


519 


Bo,„. 
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The results for the 4-31G basis set are rather similar to those of the larger basis set 
and are presented in the lower section of Table |H The combination bands are in somewhat 
better agreement than the difference bands, but overall the small basis set results are in 
good agreement with those from the larger basis set. Therefore we also computed the 
overtone, combination and difference bands in neutral coronene with the 4-3 IG basis set 
and those results are given in Table O While coronene has D6h symmetry, we use D2h 
symmetry labels to facilitate comparison with naphthalene. In this way the similarities 
between the vibrational symmetries are more readily apparent. It is important to note 
that coronene has duo hydrogens, not quartet hydrogens as in naphthalene. The CHoop 
band in coronene falls at 11.538 /im (867 cm~^) whereas in naphthalene it is at 12.677 fim 
(789 cm~^). Consequently, coronene has an overtone band in the 5 to 6 yUm region derived 
from this strong fundamental. Extrapolating to other PAHs, on the basis of typical band 
positions for CHoop bends, it is clear that PAHs containing solo and duo hydrogens will 
have overtone bands in the 5 to 6 fim region, a result consistent with the experimental 
conclusions in Sect. 13.3.11 The second point to note from the compute coronene tabulation 
is that all of the combination bands involve CHip and CHoop bending modes coupled to 
the one strong out-of-plane mode. As for naphthalene, the difference bands arise from a 
coupling of a strong CH stretch with forbidden CHip bending modes. 

The naphthalene cation was also studied using the 6-31 IG** basis set and the results 
are given in the middle section of Table HI with visualizations of the modes in Fig. [81 In 
contrast with neutral naphthalene, the first point to note is that, since ionization increases 
the intrinsic strength of many bands, the cation has several strong bands that can contribute 
to combination bands. For naphthalene, the contribution bands are mostly CHip and CHoop 
bending in character, but some begin to have contributions from modes with CC stretching 
and bending character. For example, the bands at 23 and 21 fim (435 and 476 cm~^) have 
significant CC contributions to the mode. As the cation grows in size and the number of low 



Fig. 7. — Example of one fundamental CH vibrational mode illustrative of the combination 
and difference bands of neutral naphthalene listed in Table HI The carbon skeleton is depicted 
in green, the hy ' 
by the red vect 
of the Journal. 



Fig. 8. — Example of on fundamental CH vibrational mode illustrative of the combination 
and difference bands of the naphthalene cation listed in HI The carbon skeleton is depicted in 
green, the hydrogen atoms in white and the vibrational motion of the molecule is indicated 
by the red vector arrows. Visualizations for all modes are available in the electronic edition 
of the Journal. 
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frequency modes increases, we presume that the character of the vibrations that produce 
combination bands in the cation could become even more diverse. The second difference 
between the neutral and cation forms is that the cation difference bands involve an allowed 
CH bend and forbidden CH stretches. That is, the same strong CHip bend that contributes 
to the combination bands is also the allowed band contributing to the difference bands. 

Comparing the position of the computed combination bands with the frequencies of 
the two fundamentals that comprise the bands shows only a small shift (2.6 cm~^) for the 
region of interest. This, couple d with the spec troscopic property that difference bands fall 



at (z/j — Uj) without any shift (iHerzberg 



19451 ). permits using scaled harmonic frequencies 



from larger PAHs as a guide in assessing the fundamental bands that can contribute to the 
observed overtone, difference, and combination bands. Considering the modes in three very 
large regular PAHs, C110H26) C112H26) and Ci3ofi28) we find that the cation, neutral, and 
anion forms of all three molecules should each have strong overtones near 5.4 - 5.6 fim (1818 
cm~^). Considering that we are computing the overtone as twice the scaled harmonics, the 
5.4 - 5.6 /im bands are consistent with the astronomical 5.7 fim feature as arising from 
overtones of large symmetric PAHs. All three charge states of all three molecules have 
many possible combination and difference bands in the 5-6 fim region. For the neutral 
combination bands, the strongest allowed component are CHoop bends, the same as found 
for the smaller molecules. However, it is interesting to note that there are combination 
bands in this region that are derived from an allowed CHip bend with a forbidden CHoop 
bend and an allowed CC-stretch with a forbidden drum head mode; the intensity of the 
allowed CHip bend and CC stretch is about 25% of the CHoop bend. For anions and cations, 
the intensity of the CHjp and CC stretch can exceed the strongest CHoop bending vibration. 
For large PAH cations, the CH-stretch intensity can be sizable, therefore difference bands 
in the 5-6 fim region can arise from an allowed CH stretch as well as from an allowed CH 
bend. Overall, the anion spectra in the 5-6 /im region is more similar to that of the cation 
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than that of the neutral. 

Summing up, the computational work supports the idea that PAH neutrals, cations, 
and anions have overtones arising from CHoop bends of solo and duo hydrogens that fall 
in the 5-6 yum, as do combination bands arising from in- and out-of-plane CH bends. In 
addition to the combination bands, difference bands, involving CH stretches and CH bends, 
can also fall in this region. Since there are more strong bands in the cations and anions 
than in the neutrals, more combination bands arise. It is likely that in large PAH ions, not 
all of the combination bands arise from only CHip and CCoop bends. 

For now our theoretical study is limited by the inability to compute intensities. 
Without intensities the theory is no up to making a one-to-one comparison with the 
laboratory spectra and it is impossible to definitively establish which combination and 
difference bands could contribute to the astronomical spectrum. However, the calculations 
support the idea that neutrals, cations, and anions of large PAHs can all contribute to the 
observed features. Specifically: 1) the CHoop solo and duo hydrogens lead to an overtone 
in the region of interest, 2) combination bands arising from in- and out-of-plane CH bends 
can also contribute to the observed bands, and 3) difference bands involving CH stretches 
and CH bends can also contribute to the observed bands. 



4. Astrophysical Considerations 



First we consider the 5.25 Atm band. The few ear 



its general attributes. 



AUamandola et al. 



(11989bh and 



ier studies of this 



3and focused on 



Roche et al. 



(119961 ) attributed this 



feature to overtones and combinations involving the lower lying CHoop bending modes. 
For three of the four objects considered here the 5.25 fim emission band has a FWHM of 
0.11 fim (40 cm~^). It may be slightly broader in NGC 7027 (0.14 /im, 50 cm~^), but the 
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lower signal to noise in this object makes this value less reliable. iRoche et al.l (119961 ) noted 
the striking similarity between the profiles of the 5.25 fim band and that of the 3.4, and 
11.2 fim bands. This comparison is revisited in Fig. [HI comparing the 5.25 fim band with 
the 3.4, 6.2, and 11.2 fim emission bands from NGC 7027 using the higher quality ISO 
SWS data. The major astronomical features at 6.2 and 11.2 fim arise from fundamental 
vibrations whereas the 5.25 fim band is produced by combinations and overtones involving 
fundamentals. The similarities between the profiles of these bands are striking, even down 
to the blue satellite features for the 6.2 and 11.2 fim bands. However, on an (absolute) 
energy scale the widths and skewness of the profiles show some variations, with the 6.2 /im 
profile clearly different. Both width and skewness are tied to the anharmonicity pa rameter 



19951). While 



of the transition involved, which depends on the involved mode (iJoblin et al.l 
presently relevant spectroscopic data is lacking, we note that future experimental studies 
on the relevant modes may prove to be very insightful. 



Roche et al 



(119961 ) found that the ratio of the 11.2 fim band intensity to that of the 
5.25 /xm band was roughly constant (28 - 32) over their sample of three PNe. Here, from 
Table [3] the I11.25/I5.25 in HD 44179 and the two positions in Orion are about 10 (± 1.0), 
whereas it is 14 (±1.5) in NGC 7027. We have developed a simple model to computationally 
determine this ratio based on the average (~ 9.5), experimentally established, (10 - 15) / 
5.25 yum cross sections of 15 neutral PAHs with 28 - 50 carbon atoms. The model predicts 
the I11.25/I5.25 ratio as a function of PAH size and photon excitation energies. The thermal 
approximation is us ed to determine the emission spectr um and radiative relaxation is 



taken into account (iPwek et al. 



1997 



Bakes et al. 



2OOII ) . Figure [lO] shows the predicted 



ratios as a function of PAH size between two different excitation energies, 8 and 10 eV. 
For average photon excitation energies of 8 and 10 eV, PAHs with about 40 to 65 carbon 
atoms produce the observed ratios for HD 44179 and the two positions in Orion. NGC 7027 
implies that PAHs with between 65 to about 105 carbon atoms are required to match the 



11.0 11.2 11.4 11.6 11.8 




Fig. 9. — Comparison between the normalized profile of the 5.25 /xm band to those of the 
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observed ratio. This is ver y much in hne with the estimate based on the 3.3/11.2 fim ratio 



(lAUamandola et al 



mm 



The analysis of the laboratory data and the quantum-chemical calculations point to 
the involvement of CHoop modes for the emission between 5-6 fim. The good correlation 
between the 5.25/6.2 with the 11.2/6.2 fim band strength ratios (R = 0.97) and the 5.7/6.2 
with the 11.2/6.2 fim band strength ratios (R = 0.97), which traces the CH modes, affirm 
this (Fig. [TT]) . In contrast, there is no good correlation between the 5.25/11.2 with the 
6.2/11.2 fim ratios. Showing again that the 5.25 fim and CC modes are not intimately 
connected in the astronomical spectra. 

Ionization of PAH molecules can influence band positions, but not as notably as the 
effect on the relative band intensities. Upon ionization the strength of the CC modes 
increases significantly, while the CH modes tend to decrease in strength. However, there 
appears to be a size dependence to this. PAHs of a comparable to those imp lied here do 



20081). NGC 7027 



not have significant reduction in CH band strength (iBauschlicher et al.l 
and the two positions in the Orion Bar probe the two extremes in ionizing environments. 
Future studies on a larger sample will have to establish whether the difference in 5.25/11.2 
fim ratio for NGC 7027 is a real property of the interstellar PAH spectra. While the 
theoretical data suggest that ionized PAHs and non-CHoop modes may be involved, for now 
we conclude that the correlation in Fig. [11] advocates for neutral carriers dominating the 
astronomical 5.25 fim PAH band. 



Hony et al 



(I2OOII ) made a study of the 10 - 15 fim region, the bands which are due to 
aromatic CHoop vibrational modes. Laboratory data and quantum-chemical calculations 
show that the peak position of the bands in this region depends on the number of adjacent 
hydrogen atoms on the periphery of PAH molecules. Therefore, ratios of these bands will 
reveal the structure of the emitting PAHs in astronomical environments. The predominance 
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Fig. 10. — The I11.25/I5.25 ratio as a function of neutral PAH size modelled using the thermal 
approximation and including radiative cooling. The solid line is the ratio for an average 
photo-excitation energy of 8 eV and the dashed line for an average photon energy of 10 eV. 
The black area is where the astronomically observed ratio of ~10 (±1-0) falls and the area 
in gray that where the observed ratio of ~14 (±1.5). 
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Fig. 11— Band strength ratio correlations, o HD 44179; A NGC 7027; □ Orion Bar 
H2S1; X Orion Bar D5; o Nearby galaxies (Smith, J. D., private communication). Note that 
for the Galaxies the strength of the 5.25 fim band was determined by assuming a slightly 
asymmetric P rude profile with a wid th and wavelength center trained on a high S/N subset 



of the SINGS ( iKennicutt et al. 



20031 ) spectra. Integrating the full profile produces the used 



bands strength. Also note that some of the SINGS spectra have modest redshifts, which 
bring slightly more of the feature in view, up to roughly a quarter of the feature's FWHM. 
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of the 11.2 over 1 2.7 ixm band inte nsity in the four sources investigated here points to 



compact carriers (IHony et aL 



20011 ). Since the emission in the 5-6 /im region involves 
couphng with, and intensity borrowing from the CHoop modes, this also holds for the 5.25 
and 5.7 fim bands. However, it is prudent to look at the 5-6 /im region of sources in 
which the 12.7 is as strong as the 11.2 fim band (e.g. IRAS 18317). Unfortunately good 
signal-to-noise spectra of the 5.25 /im band are not yet available for the sources known to 
show such ratios. 

Lastly, the presence of quartets produces emission in-between and beyond the positions 
of the observed astronomical 5 .25 and 5.7 Atm fe atures excludes them as important part of 



the astronomical PAH family flHony et al. 



20011 ). This is completely consistent with the 



conclusions about the PAH structures from the 10 - 15 /xm region. 

We turn now to the 5.7 fim band. Figure [6] illustrates the dependence of the 5.7 /im 
peak position on the molecular edge structure. As the laboratory data show, duos and trios 
cluster on the blue side whereas the presence of solos produce features on the red side of the 
feature. While the four spectra are quite similar, the spectrum of the Orion Bar position 
D5 seems to indicate that the duo and trio hydrogens are as important as solo hydrogens, 
even though both Orion Bar positions show similar 11.2/12.7 /im band strength ratios. 
The difference in the spectra between Orion Bar positions D5 and H2S1 suggest a different 
PAH population at both positions. This difference might be driven by erosion, where more 
irregular shaped PAHs are present at position D5. 

All of the fifteen initially considered astronomical spectra show emission features at 5.25 
and 5.7 /im, but nothing significant in-between. This rules out a large population of quartet 
containing PAHs. The spectra span a significant range of astrophysical environments, 
suggesting that many, if not most, of the astronomical PAH population is dominated by 
large and compact PAHs. It bears repeating here that the greater intrinsic strength of the 
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solo CH over that of duo CH is largely responsible for the predominance of the 11.2 /im 
astronomical band in the CHoop region. 



It has been shown by 



Peeters et al. 



(120021 ) that the variations in peak position and 



shape of the PAH bands in the 6-9 /im region can be grouped into distinctive classes. 
Moreover, these classes seem to correlate with object type. All ISM-like sources belong to 
class A. These are characterized by a CC-stretch band peaking at 6.2 fim and the '7.7' /im 
band complex peaking at 7.6 /im. On the other hand, isolated Herbig Ae/Be stars, together 
with a few Post- AGB stars and most PNe belong to class B. These are characterized by 
a CC-stretch band peaking between 6.24 and 6.28 /im and '7.7' /im band complex peaking 
between 7.8 and 8.0 /im. HD 44179, NG7027 and the two positions in the Orion Bar 
cover both of these classes, with HD 44179 and NGC 7027 representing class B and both 
positions in the Orion Bar representing class A. Interestingly enough, for the quality of 
the spectra available, the 5.25 and 5.7 /im features do not allow for such a classification, 
behavior consistent with CH modes. Given that the shift between class A and class B type 
of spect ra seems strongly influen ced by the relative populations of PAH cations to PAH 



anions (IBauschlicher et al. 



20081 ). the apparent independence of the astronomical 5.25 and 
5.7 fim bands from these classifications strongly suggests again that they track the neutral 
PAH population. 

The spectrum lab eled 'Galaxi e s' in Fig. [3] presents the slightly redshifted average 



galax y spectrum from 



(IRS; 



Houck et al. 



Smith et al. 



(l2007l ). obtained using Spitzer's Infrared Spectrograph 
2004J ) and shows part - due to IRS' cutoff at 5.2 /im - of the 5.25 /im 
band and the 5.7 /xm band. Compared to ISO, Spitzer has limited spectral resolution 
but superior sensitivity. The spectrum demonstrates the accessibility of this region using 
Spitzer and therefore the possibility for band strength analysis and, for the 5.7 fim band, 
studies of erosion processes as outlined in this paper. 
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5. Summciry and Conclusion 

This paper presents a study of the two minor PAH features in the 5-6 /xm (2000 - 1667 
cm~^) region, centered at positions near 5.25 (1905 cm~^) and 5.7 fxia (1739 cm~^). These 
contain information about the interstellar PAH population and conditions in the emission 
regions that both complement and extend the information revealed by the major bands. 

Fifteen high quality ISO SWS spectra have been investigated for emission in the 
this wavelength region, with four spectra having sufficient signal-to-noise to allow for an 
in-depth analysis. Combined with a spectral database comprised of laboratory studies 
and dedicated quantum-chemical calculations, these spectra allow us to probe the main 
characteristics of the carriers of the astronomical PAH features. After continuum and 
emission line removal, all four astronomical spectra show similar, almost universal profiles. 
However, the signal-to-noise level can be improved and there are hints of subtle, but 
interesting, variations. 

The absence of bands between 5-6 /xm in laboratory spectra of deuterated polycychc 
aromatic molecules, as well as the absence of fundamentals in the quantum-chemical 
calculations in this region along with the strong correlation between the 5.25 and 5.75 /im 
band strength with the 11.2 //m band strength in the astronomical spectra, substantiates 
the involvement of CHip and CHoop bending vibrations. In-depth analysis of the laboratory 
spectra and quantum-chemical calculations show that the astronomical 5.25 and 5.75 /im 
bands are a blend of combination, difference, and overtone bands, involving CHip and CHoop 
bending and stretching modes and, likely for the larger ionized PAHs, CCoop modes. When 
it becomes possible to compute the intensities of overtone and combination bands, this 
work should be extended. 

Turning to the hydrogen adjacency classes, PAHs with solo and duo hydrogens 
consistently produce prominent bands in the appropriate wavelength regions, whereas PAHs 
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with higher adjacency hydrogens show far richer spectra. These produce bands in-between 
and beyond the 5.25 and 5.7 /xm bands, ruhng such species out as important members 
of the emitting population. The 5.7 /im feature in itself - through its profile - contains 
adjacency class information that might be more easily accessible than through the 10 - 15 
/im CHoop region, where it is difficult to separate the duo and trio hydrogen modes. 

The data suggest that the emitting astronomical PAHs are mostly large (50 ~ A^c — 
100 ), compact, and not fully deuterated. Furthermore, both the quantum-chemical 
calculations and the absence of a correlation between the 5.25/6.2 and 5.7/6.2 fim band 
strengths with the 11.2/6.2 /im band strength ratio suggest that the 5.25 and 5.7 /im PAH 
band do not trace ionization and are carried predominately by neutrals. This point is 
reinforced by the lack of connection with class A and class B PAH band behavior. This 
suggests that high quality spectra from 5-10 /xm provide insight into the neutral as well as 
the cation and anion members of the emitting astronomical PAH family. Even so, a bigger 
collection of spectra of large compact PAHs, with mixtures of solo, duo and trio hydrogens, 
is still needed to firm up these conclusions. 
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